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The x-ray near-held speckle-scanning concept is an approach recently introduced to obtain ab¬ 
sorption, phase, and dark-held images of a sample. In this paper, we present ways of recovering 
from a sample its ultrasmall-angle x-ray scattering distribution using numerical deconvolution. We 
also show how to access the 2D phase gradient signal from random step scans, the latter having 
the potential to elude the hat-held correction error. Each feature is explained theoretically and 
demonstrated experimentally at a synchrotron x-ray facility. 


I. INTRODUCTION 

When X-ray imaging was discovered, it took only a few 
months to the scientihc community to express a great in¬ 
terest, as acknowledged by the award of the hrst Nobel 
Prize to W. Rontgen. Since then fast and tremendous 
progress could be observed in diverse applications over 
domains including medicine, material science, paleontol¬ 
ogy and many others. Within more than a century great 
minds realized the importance of X-ray imaging and con¬ 
tributed to its development with, for instance, the intro¬ 
duction of the dose concept Si and the advent of com¬ 
puted tomography ii. In parallel, constant progress in 
instrumentation allowed the development of X-ray imag¬ 
ing systems with ever better contrast, resolution and ef- 
hciency 

While questing for means of collecting more informa¬ 
tion and reducing the dose necessary to image a speci¬ 
men, the appearance of the hrst coherent X-ray sources 
favoured the emergence of phase contrast imaging meth¬ 
ods @5 S @]- Indeed, these were motivated by the fact 
that, in the X-ray domain, the material optical index 
n = 1 —5 —ij3 presents usually a refractive index 5 with a 
magnitude several orders larger than its absorption coun¬ 
terpart (3 [T^. Thus, since the inception of coherent third 
generation synchrotron sources, different near-held tech¬ 
niques of X-ray phase sensing emerged, either inspired 
by their analogue in other spectral domains or specially 
developed 0. One may categorize those imaging tech¬ 
niques depending whether they are sensitive to the Lapla- 
cian or the gradient of the phase. The propagation-based 
methods, which belong to the hrst category, make use 
of the property of contrast enhancement of edges when 
increasing the distance to the detector from the speci¬ 
men under study [iMl. The second category brings to¬ 
gether intruments such as dehection-based instruments, 
including X-ray grating interferometers [I3, Hart¬ 
mann like sensors [B Hi, coded apertures |2QL M, 
analyser based systems [22| and speckle based methods 
m, lU- Other full-held imaging approaches available 
in the X-ray domain and using different principles in¬ 
clude for instance ptychography and Zernike phase 
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microscopy [^. [^. 

Propagation based methods offer the advantage of a 
high sensitivity to sharp changes in electronic density, 
the sample high frequency features being rendered with 
high contrast thanks to the phenomenon of optical inter¬ 
ference. Conversely, they often suffer from reconstruc¬ 
tion artefacts when imaging materials presenting a slowly 
varying density and low spatial frequencies due to the 
very nature of the recorded signal, i.e. a Fresnel diffrac¬ 
tion image [28|. In practice, a phase gradient sensitive 
method would be preferable to image a homogenous ma¬ 
terial, although usually resulting in a spatial resolution 
lower than the one obtained with a propagation based 
technique, all other conditions remaining equal. For 
their implementation, phase-gradient sensitive methods 
require the use of a specific optics (grid, grating, crys¬ 
tal...) to modulate the wavefront, which can, sometimes, 
be expensive or difficult to get. 

Darkfield X-ray contrast is an additional contrast mode 
inspired from its counterpart in the visible light and made 

a sible thanks to the availability of crystal analysers 
[ 29 I and of grating interferometers Q. Darkfield 
imaging permits to reveal the structure and the orienta¬ 
tion of features having a characteristic size much smaller 
than the detector resolution. For instance, the map¬ 
ping of the dark field signal with a grating interferom¬ 
eter allows one to determine the local reduction of the 
interference amplitude caused by X-ray beam propaga¬ 
tion through the sample. The interference amplitude ac¬ 
tually reflects the progressive degradation of the beam 
transverse coherence properties caused by scattering in¬ 
duced upon photon propagation through the samples 
sub-microscopic features. 

The X-ray near field speckle scanning technique [3l|,[32| 
is a recent and particularly attractive deflection based 
technique as it provides both phase and darkfield con¬ 
trast images in addition to the traditional absorption 
one. Moreover, it uses only a random phase object as 
wavefront modulator. This phase object is used to gen¬ 
erate a speckle intensity pattern in the near field domain 
caused by mutual interference of the scattered waves. 
The near-held domain is particularly interesting because 
the speckle pattern distortion is only ruled by the wave- 
front transformation upon propagation in the Fresnel 
zone [ 33 I, [U , a region easily accessible with X-rays [35| . 
The X-ray near held speckle tracking technique pre- 
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FIG. 1. (Color online) First image recorded during the scan 
with (a) and without an object (b) inserted in the beam, (c) 
Experimental setup for the data collection. 


viously developed and easily implemented mil could 
provide phase information from one single exposure of 
the sample. However, its resolution is limited by the 
speckle grain size and artefacts may appear in the case 
of complex objects. The speckle scanning method ap¬ 
peared to overcome these limitations, nonetheless at the 
cost of additional image acquisitions. The technique ad¬ 
vantages include a high sensitivity to the phase gradient, 
a moderate requirement for the beam monochromatic¬ 
ity, a high photon efficiency (no absorption in the sys¬ 
tem) and the use of a minimalist optical element in the 
form of a scattering membrane. Albeit the speckle based 
methods were first developed at a synchrotron source, 
they were recently adapted with success to a laboratory 
source, thereby broadening the range of applications . 
In this letter, several new aspects of the speckle scanning 
concept are presented. After having recalled the basic 
concept and the general experimental setup, we will ex¬ 
pound a method capable of accessing the ultra-small an¬ 
gle X-ray scattering distribution of a sample from 2D 
rasters scans. Next, a processing scheme to recover the 
phase gradient information in two dimensions from scans 
with random steps is demonstrated. Each aspect is illus¬ 
trated with experimental applications. We also describe 
a way of eluding the ’fiat-field correction error’ using the 
speckle scanning technique. A short discussion dealing 
with the sensitivity and operational constraints is con¬ 
cluding this paper. 


II. EXPERIMENTAL SETUP 

Let us introduce the experimental setup that was used 
to acquire the data presented further. The experiments 
were conducted at the beamline BM05 of the ESRF • 
There, the photons are produced by synchrotron radia¬ 
tion from a bending magnet of 0.85 T installed on a stor¬ 
age ring operating with 6.04 GeV electrons. In the exper¬ 
iments presented, a beam with an energy E = 17 keV was 
selected thanks to a double crystal Si(lll) monochroma¬ 
tor with a selectivity of AEjE ^ 10“^. The experi¬ 
mental setup sketched in Fig. [T] was installed in a lead- 
shielded hutch, at distance from the scattering object to 
the source of R ^ 55 m. Considering the source size of 
29 jiim rms vertically and ^77 jam rms horizontally, the 
transverse coherence at the level of the membrane at E 
= 17 keV is expected to be at most ^29 /im vertically 
by 11 jam horizontally. 

The scattering object, a piece of sand paper with a 
grit designation P800, was mounted on a motorized 2D 
piezo stage, itself mounted above stepping micromotors. 
Such a setup permitted to move the random phase object 
transversally to the beam with either a nanometer preci¬ 
sion and a 100 jam range, or a micrometer accuracy over a 
millimeter range. The scatterer was placed at a distance 
I = 450 mm upstream of the sample and the detector at 
a variable distance d downstream from it. The imaging 
system was a FReLoN (Fast Read-Out Low Noise) CCD 
camera receiving a visible light signal resulting from the 
luminescence produced by X-ray illumination of a scintil¬ 
lator and imaged by a microscope objective. The system 
effective pixel size was of Sp = 5.S jam. Near field speckle 
of a size of a few pixels was observable in the recorded 
images with a contrast defined by a coefficient of varia¬ 
tion of ^ 14% when d = 850 mm (visible in Fig. [1] (a) 
and (b)). 


III. SCANNING TECHNIQUE WITH 
CONSTANT STEP 

This speckle scanning mode is mainly presented in 
Ref. [3l| where the conceffi is described as the crossroads 
of the speckle tracking [^, and grating interferome¬ 
try techniques [^, . Since we are going to present an 

extension of it, let us recall quickly the method principle. 

Figure [2] shows the geometry considered wvith the vec¬ 
tors (ex,ey) perpendicular to the beam propagation di¬ 
rection Gz. The scattering object with random phasors 
is placed in a partially coherent X-ray beam, either right 
upstream or downstream the investigated object. While 
small differences exist between these two configurations 
1110 they are negligible in this following experiments 
as the beam was nearly collimated. Nonetheless, we will 
rederive here the correct formulae for the magnifying ge¬ 
ometries currently encountered, e.g. when using a focus¬ 
ing optics or a highly divergent laboratory X-ray source. 
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FIG. 2. (Color online) (a) Geometry for the scanning X-ray 
speckle imaging technique with constant steps, (b) Example 
of a reference speckle pattern built from the intensity values 
read in the yellow pixel during the scan, (c) Corresponding 
pattern recorded in the same pixel during a scan with the 
sample in the beam, (d) Scattering distribution obtained by 
deconvolution of (b) with (c). 


The method consists of coupling together data col¬ 
lected when raster scanning the membrane with the sam¬ 
ple present in the X-ray beam (e.g. array of Fig.[Tl(c)) or 
removed from it (Fig. [11(b)). For the technique to be ef¬ 
fective, the membrane must be scanned with a positional 
repeatability p = {X,Y) of a fraction of a speckle grain 
size and with a strictly constant step size Ss = cst. A 
pair of patterns recorded in a pixel located at r = (x^y) 
is considered to be independent from the data collected 
in the neighboring pixel. The pixel marked in yellow in 
Fig. ISfa) with the pattern pair displayed in (b-c) is an 
example of data sets recorded at a given position r. Then 
the 2D intensity value array recorded when scanning the 
membrane with the sample in the beam path and noted / 
(Fig.|2l(c)) is cross-correlated with the one obtained dur¬ 
ing the reference scan and noted g (Fig. |2j(b)), i.e. with 
no sample in the beam. The outcome of the operation is a 
correlation map whose maximum peak location indicates 
the displacement vector v between the two arrays: 

v(r) = argmax / /(r,p)5(r,P + r)dr (1) 

Well known recipes based on correlation coefficient 
curve-fitting permit to calculate v with a fraction of a 


step accuracy [4l|, [42[. Then, from the angles a, 7 and ( 
defined in Fig. HJand for the two geometries (a) or (b), 
we have a = 7 + and so by calculations: 


a 


RYd^l 

R 



( 2 ) 


which is valid under the small angle approximation en¬ 
countered with X-rays. For a collimated beam, F = 
{Rd 1)/R reduces to 1 and (a = C- 
In the following we use the del operator noted V. As 
the angle a relates to the differential wavefront grad ient 
a = VW itself linked to the phase gradient [43, we 
can write: 


k\/W = V0(r) 


/c-^Fv(r) 


(3) 


where k is the wavenumber equal to 27r/A, with A the 
photon wavelength. Most often, data recorded through 
2D scans is followed by projection of v in order to ex¬ 
tract the two orthogonal transverse phase gradient com¬ 
ponents. One-dimensional scans with much less points 
have also proved to be efficient for recovering the ID 
phase gradient [3l|. The absorption, or here the trans¬ 
mission image T, can be correctly calculated free from 
the fiat field correction errors that occur in presence of 
a strongly dephasing object and/or of a high magnifica¬ 
tion projection geometry [25|, Ej- Its formulae can be 
obtained by consideration of Fig. jH 


T(r) = 




g{g{TYTssv)) 


(4) 


where g denotes the mean operator over the variable p. 
Finally, the darkfield image Df can be reconstructed by 
considering the decrease of speckle contrast that occurs 
when introducing the sample in the beam. More pre¬ 
cisely one has to calculate the ratio of the coefficient of 
variations of / and g\ 


^ 1 ^(/W) .5. 

r(r) C7(5r(r+ rssv)) 

with (7 the standard deviation operator over p. 

Using a blueberry as a sample, the Figure |3l( a-d) il¬ 
lustrates the various modes the technique offers. The set 
of two scans was recorded with d = 850 mm, Sg = 3 /im, 
each scan being composed of an array of 25 x 25 points. 
The collection of such large data scans was performed 
in an attempt to optimize the calculation of the higher 
moments of the scattering distribution that require more 
statistics, as we shall see further. Otherwise many less 
points would be needed. Figure [3] presents the calculated 
transmission image in (a), in (b) and (c) the respectively 
horizontal and vertical wavefront gradients VlU.ex/y, in 
(d) the integrated phase image, and in (e) the darkfield 
image. 


























4 



FIG. 3. (Color online) Image modes of a blueberry sample obtained with the method of Sec. lIIIl (a) transmission image T, (b-c) 
horizontal and vertical wavefront gradients VVF.ex/y, (d) phase reconstruction 0 and (e) scattering image Df. Images obtained 
using the method of Sec. IIVI (f-g) horizontal and vertical wavefront gradients VVF.ex/y. Higher moments of the scattering 
distribution: (i) \/M^, (k) (1) ^Mqs and (h) \^M^. 


IV. ULTRA SMALL ANGLE X-RAY 
SCATTERING DISTRIBUTIONS 

Recent works demonstrated a way of recovering the 
scattering distribution of a sample in each detector pixel 
from raster scans. The co ncep t was first introduced by 
Yashiro et al. in [45| and [46| using a grating interfer¬ 
ometer. Wherein, the scattering effect of the sample is 
treated as a local point spread function whose estimated 
width mirrors the local scattering strength. Later on, 
Modregger et al [13, |48| also used a grating interferom¬ 
eter but, this time, combined with numerical deconvo¬ 
lution to be able to extract the higher moments of the 
scattering distribution. In this section, we introduce by 
equivalence the recovery of the ultra-small angle X-ray 
scattering distribution using near field speckle scans sim¬ 
ilar to the ones recorded and used in the previous section. 

Let us consider the light distribution / in the case of a 
sample present in the X-ray beam path. Upon propaga¬ 
tion, and for a given pixel position r, we can consider / as 


the effect of the sample optical transfer function o(r, p) 
convoluted with the input reference signal g [^, : 

/(r, p) = y p “ ^) 9 {r, p)du (6) 

Henceforth, for each pair of patterns (with and with¬ 
out sample) (/, 5 f)(r), a deconvolution procedure can be 
applied to retrieve the distribution o. For this p ur p ose 
the iterative Richardson-Lucy deconvolution |49l-[^ al¬ 
gorithm is convenient because it is robust and its input 
does not require any assumption on the system noise . 
Such processing was applied to the blueberry sample data 
of Sec. mil with 100 iterations used in the Richardson- 
Lucy deconvolution for each pixel. An illustration of the 
scattering distribution o resulting from the deconvolution 
of the data of Fig. [2j(b) and (c) is displayed in (d). 

Geometrically, the position p corresponds to the an¬ 
gular coordinate a = F^p (see Fig. m). By projection 
in this angular variable coordinate system, a scattering 
image o{a) can be constructed for each angular vector. 
























































5 


(a) 

Membrane 

Detector 

_ 

-— 

Source 



◄- 

m ScHUplc 

-- 


(b) Detector 


Sampl 

Membrane 

e ; . . 


_ ^ 

-r . . T ^ 

|v r . . 


Source 



1 

◄-► 


-► 


FIG. 4. (Color online) The two possible geometries for the 
technique with the membrane placed either (a) upstream or 
(b) downstream the investigated object. The refraction angle 
of the (dashed) red ray passing through the sample is noted 
Oi and the (plain) blue ray shows the ray falling in the same 
detector pixel when no sample is in the beam. 




FIG. 5. Array of the 5 X 5 central scatter images organized 
by their angular vector a — Ssp/d . 


Figure [5] shows an image array of the blueberry sample 
organized as a function of the scattering angular coordi¬ 
nates. Given a scan step of 5 ^ = 3 /im, the angular step 
size from image to image was of 3.5 /irad. 

Moments of the scattering distribution contain phys¬ 
ical information about the sample features and corre¬ 
spond, as we shall see for the first of them, to the tra¬ 
ditional contrast modalities, i.e. the absorption, phase 
and scattering signals. The central moments M^an of or¬ 
der + n > 1 of o(r) are calculated within the 


meaning of angular probability density functions: 

1 / s \ r _ _ 

Mmn{r) = — (^) j {X-X,r(Y-Yrroir, p)dp 

(7) 

where: 

Moo = J o{r, p)dp (8) 

and (Xr,Fr) are the centroid positions of o(r, p): 

— _ JXo{r,p)dp - _ fYo{r,p)dp 

. , hr — (y) 

Moo 

We have then v(r) = (Xr, by definition of v, from 
which we can calculate using Eq. [3l Furthermore, 
the two second normalized moments \/Mq 2 and 
mirror the angular characteristic scattering width in the 
two transverse directions. Higher moments provide com¬ 
plementary information on the scattering distribution, 
with, for instance, the third and fourth moments related 
to its skewness and kurtosis. The wavefront differential 
gradients VW.e^/y for the blueberry sample obtained 
using the scattering distribution centroid are shown in 
Fig- [31(f) and (g). Even visually they are in good agree¬ 
ment with the equivalent maps (b) and (c) calculated 
with the constant step method presented in the previous 
section. In this same figure, the second moments y/M^ 
and \/M 2 o are shown in (i) and (j) and the higher mo¬ 
ments ^Mso, ^Mo 3 and drawn in (k,l) and (h). 

V. SCANS WITH RANDOM STEPS 

In this section, we present a processing method to re¬ 
cover the sample-induced 2D phase gradients from im¬ 
ages acquired using random-step scans, i.e. scans where 
the images are recorded with the scattering membrane 
moved at random transverse positions (X, Y) in the plan 
(gx, Gy), with no particular translation distance from one 
image to the next. However, it is crucial that scans per¬ 
formed with and without the sample have their images 
taken with the scattering membrane located at the very 
same positions for each corresponding step. 

Here, we will change slightly the notation, in order 
not to confuse the reader with the previous process¬ 
ing method. We note i(r, n) the intensity collected at 
the position r of the image of a scan containing 
N images when the sample is inserted into the beam. 
Then, writing ir,n = n), we have the vectors ir = 

(ir,i, - •••5 '^r,N) built from N realizations of the ran¬ 

dom variables Jr, induced by the variable position of the 
membrane. With equivalent notations, we also build the 
reference vectors oor, from N realizations of H4 obtained 
in absence of sample. 

The idea is here again to use correlation calculations 
between the speckle distribution in the various pixels to 
retrieve the beam wavefront derivative through the calcu¬ 
lation of the local light deflection angles (a(r) as sketched 
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FIG. 6. (Color online) Geometry of the concept, (a) Data 
vector ir built from N images with the sample in the beam, 
(b) Set of reference vectors {ct;r+h} in the reference data, (c) 
Correlation factor map of p(ir, Ct;r+h)- 


in Fig. [6l So, we use the Pearson correlation coefficient 
Pr generally defined by: 


p{Ir,Wr,) 


E[{Ir-p{Ir)){Wr,-p{Wrr))] 

a{Ir)a{Wr,) 


( 10 ) 


where E denotes the expected value. The factor p is the 
normalized covariance of the pixel signal distributions. It 
is used as figure of merit to evaluate the level of similarity 
between the data acquired in the different pixels of the 
two scans. For the vectors i and cj, this operation can 
also be seen as a scalar product. Thus, we calculate an 
estimate of each projection vector of ir onto vecrors cjr+h 
located at a small distance h from r in the reference 
scans: 


N _ 

{E,k '^r)(^r+h,/c ^r+h) 

p(ir,CJr+h) = / 

N _ ^ 

4'^r)^ (^r+h,/c ^r+h)^ 

V /c k 

_ ( 11 ) 
where E and oor denote the mean values of the measure¬ 
ment vectors. 

For all ir, and with the correlation map p(ir, c^r+h) (see 
Fig. [61(c)), we can locate the vector position of c^r+h^ax 
rendering maximum correlation: 

h-maai = argmaxp(ir,Wr+h) (12) 

h 

Here also the vector h^ax is determined with a sub¬ 
pixel accuracy by numerical interpolation over the neigh¬ 
boring pixels [^, . Following the law of light propa¬ 

gation in homogeneous media, the phase gradient V(/) is 


recovered by: 


V<^(r) = (13) 

where Sp is the detector pixel size. The corrected trans¬ 
mission signal is equal to: 


Pjlr) 

P{Wr+h^^^(r)) 


(14) 


and is valid for any geometry. With this processing 
method the darkfield signal writes: 


Df 


1 <T{Ir) 

Tr (^(W^r+h„a.(r)) 


(15) 


which is fully equivalent to Eq. |5l This method can be 
seen as an over-sampling version of the X-ray speckle 
tracking technique. Indeed, the cross-correlation opera¬ 
tion is operated between data collected in different pixels 
as performed with speckle tracking [^. Moreover, the 
refraction angle a is here calculated at the sample posi¬ 
tion (see Fig.|6j) and not at the detector plane level as in 
the case of the method of Sec. uni with the angle (. 

An experimental demonstration of the processing 
scheme is provided in Fig. [71 using the same blueberry 
sample as before, for comparison purpose. A pair of 25- 
image scans was acquired by applying displacements to 
the membrane with steps larger than 20 pm using a mi¬ 
crometer precision linear positioning motor. Fig. [71(a) 
and (b) are respectively showing the horizontal and ver¬ 
tical transverse differential wavefront gradients. The re¬ 
sults obtained are in nearly perfect agreement with the 
one calculated with the constant step size method pre¬ 
sented in Sec. IIIIl Figure [71 shows the reconstructed phase 
shift of the sample (c), the transmission image (d) and 
the darkfield image (e). Figure [71 (f) displays a plot of 
cuts along the marked lines of the three previous images 
normalized to an arbitrary scale. One can observe the 
way the different contrasts vary depending on the image 
modality. 

This processing scheme is particularly interesting as it 
allows to recover the 2D phase and scattering informa¬ 
tion for every pixel, as does the scheme with constant 
step, although requiring a lower number of images. With 
the constant-step method, 0{N‘^) images are necessary 
to build the 4D matrix /(r, p) whilst only 0{N) are re¬ 
quired in this method for 2D phase sensing. 


VI. DISCUSSION 
A. Sensitivity 

The sensitivity of the speckle scanning technique differs 
depending on the processing method. When applying the 
algorithm presented in Sec. [Till on data acquired with a 
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FIG. 7. (Color online) Random step processing scheme, (a) 
Horizontal and (b) vertical differential wavefront gradients, 
(c) Phase (d) transmission and (e) scattering reconstruction 
images, (f) Intensity plot along the lines marked in (c-e). 


scan of constant step, the accuracy a (a) on the wavefront 
gradient is: 

o-i(a) = ^r(T(Vr) ( 16 ) 

This can be compared to the processing method intro¬ 
duced in Sec. m for instance when using the same data: 

0'2(a) = (17) 


From these equations, with a similar experimental 
setup, it is the constant-step scanning method that de¬ 
livers the highest sensitivity. As a matter of fact, consid¬ 
ering an equivalent sub-pixel/step accuracy for the peak 
finder algorithm in both processing methods, we have 
cr(vr) ~ cr(hr) and we can set Ss < Sp to diminish ai. 
When working with a collimated X-ray beam, the gain 
in sensitivity gets rapidly limited by the minimum step 
length necessary to get a sufficient intensity variation 
within a pixel from one step to the next [^, [bJl- On 
the contrary, with a magnification projection geometry 
setup, the gain can be tremendous since a tiny displace¬ 
ment of the membrane located near the X-ray beam fo¬ 
cal point can generate a much larger displacement of the 


speckle pattern at the detector level. In this case, when 
generating a sufficient speckle statistic variation for the 
correlation algorithm to be efficient and accurate, the 
angular resolution may reach the nanoradian scale [55| . 

The convergence of the distribution estimators in 
Eq. [11] scales with 1 /\/N [56| . This implies that only 
the few first images are necessary to reach a suitable ac¬ 
curacy since ^/Na{h) 0 and \/]Vcr(v) 0 as N ^ oo. 

From Eq. [TOl considering two identical scans or areas of 
r where the X-ray light is neither refracted nor scattered 
by any sample, the distributions (Jr, Wr) should, by def¬ 
inition, render a theoretical value p(/r, Wr) = 1. Eor our 
experimental data, when using the method of Sec. El the 
deviation from unity of p over an area of 150 x 150 pix¬ 
els amounts to (j{a) = 0.8 /irad for a scan of 25 images. 
When more than 150 images are used it tends to a limit 
of (j{a) = 0.35 /irad. This last figure is comparable to 
the one calculated with the two constant-step methods 
presented in Sec. lIIIl and lIVi suggesting that the discrep¬ 
ancy may be due to monochromattor instabilities, more 
felt in the vertical direction. 


B. Tracking orientation 


Section O explains and illustrates the way the refrac¬ 
tion angle is recovered by examining the vector positions 
in the original reference scan data. This subsection aims 
at illustrating the importance of operating along this line 
rather than tracking the reference scan vectors in the 
data array collected when the sample is present in the 
beam. The problem raised here is due to the creation 
of optical vortices in the propagated beam when going 
through the sample discontinuities • 

EigurelHIpresents an experimental data set for a juniper 
berry sample imaged with X-rays. That sample exhibits 
strong scattering features and turbid phase shift. Eig- 
ure IHl(a) and (b) are successfully displaying the trans¬ 
mission and phase contrast images, respectively. Eig- 
ure IHl(c) presents the horizontal speckle displacement 
vectors calculated picking vectors within the sample scan 
and searching for their counterpart vectors in the refer¬ 
ence scan. Conversely, (d) corresponds to the tracking 
of vectors from the reference scan across the sample scan 
data. Eigure|8l(e-f) are zooming in the square subsets of 
(c) and (d), while (g) shows the absolute intensity differ¬ 
ence between these two subsets. The significant differ¬ 
ence observed between the two calculated images results 
from the incapacity of the algorithm to seize reference 
vectors folded up in optical vortices upon propagation. 
As these vectors add up at singularity points they gen¬ 
erate new distributions A [13 that are no longer corre¬ 
lated to any reference vector. Obviously, as the detector 
pixel size and the speckle grains themselves are getting 
smaller, or as the propagation distance increases, the ef¬ 
fect becomes more pronounced. 
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FIG. 8. (Color online) (a-b) Correct transmission and phase 
images of a juniper cone, (c) Horizontal phase gradient calcu¬ 
lated searching the location of vectors from the sample scan 
across the reference scan, (d) Output of the reciprocal calcu¬ 
lation, searching reference vector locations across the sample 
scan, (e-f) Zoom in the small subset drawn on the above 
images, (f) Absolute difference between the values of two 
magnified image subsets: A = (e) — (/). 

VII. CONCLUSION 

Several advanced processing schemes were demon¬ 
strated within the context of a new X-ray speckle scan¬ 


ning imaging technique. Starting from the presently 
available scanning method with a constant step, we pro¬ 
posed ways of building error free transmission images and 
of accessing the sample scattering distribution, pixel by 
pixel, through a numerical deconvolution of the scans 
data. The various processing schemes were demonstrated 
with experimental data, assessing thereby the excellence 
of the method. The recovery of the sample scattering dis¬ 
tribution is of interest in various fields of material science 
as well as in biological imaging where the orientations 
of sub-pixel structures are responsible for the functional 
behavior of the systems. In addition, a new processing 
scheme with relaxed constraint on the step size was pro¬ 
posed. It can be seen as an over sampled version of the 
speckle tracking technique with a resolution enhanced up 
to the one defined by the detector. The speckle-scanning 
imaging technique based on a random step scheme re¬ 
duces dramatically the number of acquisitions as com¬ 
pared to the previously available schemes. This point is 
often essential in the case of biological imaging, where the 
photon efficiency permits to reduce the dose delivered to 
the sample. Future work will permit to push and test 
the limit of the method in this regard. The proposed ap¬ 
proach will also promote the development of the speckle 
imaging technique with laboratory low brilliance sources. 
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